Peeling of Dirac and Maxwell fields on a 
Schwarzschild background 
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Abstract 

We study the peeling of Dirac and Maxwell fields on a Schwarzschild background following 
, the approach developed by the authors in [16] for the wave equation. The method combines 

£NJ ' a conformal compactification with vector field techniques in order to work out the optimal 

| , space of initial data for a given transverse regularity of the rescaled field across null infinity. 

The results show that analogous decay and regularity assumptions in Minkowski and in 
Schwarzschild produce the same regularity across null infinity. The results are valid also for 
the classes of asymptotically simple spacetimes constructed by Corvino-Schoen / Chrusciel- 
(N ; Delay. 
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1 Introduction 

Zero rest-mass fields on asymptotically flat spacetimes admit a peeling-off, or peeling, property, 
discovered by Sachs in 1961 [23] in the flat case. It can be described in terms of principle null 
directions : for a zero rest-mass field 4>AB...F = <fr(AB...F) with n indices, the part of the field 
falling-off like r~ fc_1 , < k < n, along outgoing null geodesies, has n — k of its principle null 
directions aligned along the generator of the geodesies. On Minkowski spacetime, this can be 
stated simply in terms of components in a well-chosen spin-frame. Consider the Newman-Penrose 
tetrad 

1 1 1 i 

I = -7=(d t + d r ), n= -7=(d t -d r ), m = ~^={d e + -r-^) , rh 
V2 V2 r\/2 sin ft 

and the associated spin-frame {o , l} (unique modulo overall sign), then the component (f) n —k^ 
which is the contraction of <p with n-h's and k o's, falls off like r~ k ~ l along the integral curves 
of I. 

It is essential to note that the property as we have brutally stated it is wrong. One can 
consider initial data on a spacelike slice that are exponentially increasing at infinity and the 
associated solution will not fall-off at all along outgoing null geodesies. The peeling is true in 
flat space-time for smooth compactly supported initial data and for certain classes of data that 
satisfy adequate regularity and fall-off assumptions. The situation is expected to be similar on 
generic asymptotically flat spacetimes but it has been speculated that the conditions on the 
initial data may be more stringent due to the more complicated asymptotic structure. 

In 1965, Penrose [21] presented a new derivation of the peeling based on conformal compacti- 
fications. Using the conformal embedding of Minkowski spacetime into the Einstein cylinder and 
the conformal invariance of zero rest-mass field equations, he showed that the peeling property 
is equivalent to the continuity of the rescaled field at null infinity (j^). Then he argued that the 
peeling should be a generic behaviour of zero rest-mass fields on asymptotically flat spacetimes. 
He went on to define a class of spacetimes, referred to as asymptotically simple spacetimes, pro- 
viding a generic model of asymptotic flatness as well as a framework in which the peeling should 
occur under reasonable conditions. These spacetimes are classifed according to the regularity of 
their conformal metric at null infinity, which encodes the information of the peeling of the Weyl 
tensor. 

His results raised criticisms regarding the genericity of asymptotically simple spacetimes and 
the peeling property. The jist of the arguments put forward was that the asymptotic structure 
of the Schwarzschild metric would only allow peeling for a considerably smaller family of initial 
data than in Minkowski spacetime. Indeed, asymptotically simple spacetimes are meant to be 
physically reasonable and as such contain mass-energy, which means that the physical metric 
differs from the flat one at first approximation by a Schwarzschild- type behaviour in m/r. The 
asymptotic structure of Schwarzschild's spacetime is substantially different from Minkowski's in 
particular with a singular conformal structure at spatial infinity. It was unclear as to whether 
this would impose more stringent hypotheses on initial data to ensure peeling; the singularity 
at spatial infinity could interact with the tail of the falloff of the initial data at spatial infinity 
and prevent the peeling that might otherwise take place if the data were compactly suppported. 
The genericity of the peeling for zero rest-mass fields being questioned, the peeling for linearized 
gravity on Schwarzschild's spacetime was also doubted and this made the asymptotically simple 
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spacetime model appear as anything but generic. The question of regularity of null infinity 
and asymptotic simplicity has now been resolved in various ways, Christodoulou-Klainerman 
[2], Corvino [5], Chrusciel and Delay [3, 4], Corvino-Schoen [6], Friedrich (see [10] for a survey 
of his contributions) and Klainerman-Nicolo [12, 13, 14]. However, even in the simple case of 
the Schwarzschild metric, it was not at all clear, until the authors provided a first element of 
answer in [16], whether zero rest-mass fields admit peeling properties for reasonably large classes 
of initial data. 

Penrose's constructions in [21] provide powerful techniques to analyze the constraints on 
the initial data implicit in the peeling property. A t = constant slice of Minkowski spacetime 
corresponds to a 3-sphere with a point removed on the Einstein cylinder. Physical initial data that 
once rescaled extend as smooth functions on the whole 3-sphere give rise, as ensured by Leray's 
theorem, to a rescaled solution that is smooth on the whole Einstein cylinder and consequently 
the physical solution satisfies the peeling property. This is a little crude since the peeling really 
means continuity of the rescaled solution at null infinity and the class of data considered above 
provides solutions that are C°° across ,f. A more detailed understanding requires to consider 
intermediate regularities, but C k spaces are not adapted to the Cauchy problem for hyperbolic 
equations ; a typical example is the wave equation on Minkowski spacetime : for initial data 
f\t=o e C k (R 3 ) and d t f\ t=0 € C fc " 1 (M 3 ), the solution / is generally not in C k (R x M 3 ). It is 
better to use Sobolev spaces instead since they are naturally controlled by energy estimates. In 
a previous paper [16], the authors used a combination of Penrose's conformal techniques and 
geometric energy estimates (or vector field methods) to characterize completely the spaces of 
initial data ensuring peeling for the scalar wave equation on the Schwarzschild metric. A new 
definition of peeling with specifyable order of regularity was given in terms of (weighted) Sobolev 
spaces. The classes of data were shown to have analogous fall-off properties as those obtained 
in the flat case using the conformal embedding in the Einstein cylinder. This established that 
at least for the wave equation, the different asymptotic structure of the Schwarzschild spacetime 
does not change the classes of data ensuring peeling. To gain a more complete understanding of 
the question, it is necessary to study other types of fields on a Schwarzschild background, such 
as higher spin zero rest-mass fields, or solutions of non linear equations, and then to extend the 
results to more general asymptotically flat spacetimes. 

The present paper uses similar techniques to investigate the peeling for Dirac and Maxwell 
fields on the Schwarzschild spacetime. It is organized as follows. Section 2 describes the geo- 
metric ingredients of the method : the conformal (partial) compactification of Schwarzschild's 
spacetime with choices of Newman-Penrose tetrads on the physical and rescaled spacetimes, the 
corresponding rescaling of spin-coefficients, the neighbourhood of spacelike infinity in which we 
establish our estimates, a choice of foliation and the Morawetz vector field. This vector field is 
as crucial for the Maxwell case as it was for the wave equation, but it does not play any part 
in the construction for Dirac fields thanks to the existence of a conserved current independent 
of a choice of observer. In section 3, we focus on the field equations, their conformal invariance, 
the proof of the equivalence between the peeling and the continuity of the rescaled field at null 
infinity as Penrose gave it in [21], the conserved quantities (current for Dirac and stress-energy 
tensor for Maxwell), the associated energies on the hypersurfaces we work with and the generic 
method for obtaining energy estimates for a perturbed equation (explained simply in the Dirac 
case). Section 4 contains the peeling results and their proof for Dirac and Maxwell fields. We 
establish energy estimates for fields supported away from spacelike infinity and then use these to 
construct our spaces of initial data by completion in the norms obtained on the initial hypersur- 
face. We get a full set of function spaces with all degrees of regularity ; for data in these spaces, 
we have estimates both ways between the norm on the initial hypersurface and a corresponding 
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norm on involving transverse derivatives as well as angular ones. The main difficulty is of 
course the estimates on transverse derivatives. The control of angular derivatives is straightfor- 
ward thanks to the spherical symmetry ; in the case of Dirac fields, it can be obtained elegantly 
by commuting into the equation the Dirac operator on the sphere ; for the Maxwell system, this 
cannot be done because we have four equations and three unknowns (this is also manifest in 
the spin and boost weights of the components of the field and their transformation under the 
Geroch-Held- Penrose angular operators), so we use a set of three Killing vectors on the sphere 
instead. The main results, given in theorems 1 and 2, are the energy estimates and the func- 
tion spaces on the initial data surface that one can infer from them. The interpretation of the 
results, namely the comparison between our classes of initial data and the classes obtained in 
Minkowski spacetime using the full conformal embedding in the Einstein cylinder, is given in 
section 5, together with a remark on the constraints. All the calculations in this paper are done 
using partial derivatives of the field components. One may prefer using covariant derivatives of 
the fields and then taking components. This requires to know the different curvature spinors of 
the rescaled Schwarzschild spacetime. The calculation of these quantities and the derivation of 
an energy estimate for a transverse derivative using this approach is given for Dirac fields in the 
appendix. Although the calculations are more involved, the error terms in the conservation law 
for the transverse derivative are much simpler than if we use partial derivatives (see remark A.l 
at the end of the paper) . 

Since the main nontrivial results of this paper concern the behaviour of fields in a neigh- 
bourhood of space-like infinity and its intersection with null infiinity, the results are valid for the 
classes of spacetimes of Corvino-Schoen / Chrusciel-Delay. 

Notations. Throughout the paper, we use the formalisms of abstract indices, 2-component 
spinors, Newman-Penrose and Geroch-Held-Penrose. 

2 Geometric setting 

2.1 Rescaled Schwarzschild spacetime 

We work on the Schwarzschild metric 

g = F{r)dt 2 - F(r) _1 dr 2 - r 2 du 2 , m > , 
F(r) = 1 - 2m/r , duo 2 = d9 2 + sin 2 Odip 2 , 

in the region outside the black-hole M. t x]2m, +oo[ r xS 2 . Introducing the variables u = t — r*, 
where r* = r + 2mlog(r — 2m) is the Regge- Wheeler coordinate, and R = 1/r, we obtain the 
following expression for the metric g conformally rescaled using the conformal factor R : 

g = R 2 g = R 2 Fdu 2 - 2dudR - du 2 , still denoting F = F{r) = 1 - 2mR . (1) 

These choices of conformal rescaling and variables allow to define naturally future null infinity 
(J f+ ) as R u x {R = 0} x S 2 . The Levi-Civita symbols must be rescaled accordingly : 

SAB = R^AB ■ (2) 

We make on the exterior of the black hole the following choice of unitary (for g) Newman-Penrose 
tetrad : 
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with corresponding dual tetrad 

[~F 1 2 I R?F \ 1 

l a dx a = y — du , ri a dx a = y — ( — ^— du - dRj , m a dx a = - -j= (d9 + i sin 9d(p) . (4) 

This is in fact a simple rescaling of a classic Newman-Penrose tetrad for the metric g : 

l a da = -^(d t + d r J, n°d a = -J=(d t -d r J, m *d a = -±=(do + ^t 

since we have 

l a = r 2 l a , h a = n a , m a = rm a . 

The indices for the rescaled frame vectors are lowered with the rescaled metric, so we have the 
following link with the unrescaled frame co-vectors : 

L = l a , n a = R 2 n a , rh a = Rm a . 

In terms of associated spin- frames, this corresponds to the rescaling 

6 A = ro A , i A = i A , 6a = oa , la = Rla ■ (6) 

We shall use the standard Newman-Penrose notations D, D' , 5 and 5' for the directional deriva- 
tives l a V a , n a V a , m a V a and rh a V a ; similarly, we denote by D, D' , 5 and 5' the directional 
derivatives along I, n, rh and in. 

The 4- volume measure associated with the metric g is given by 

dVol 4 = il An Am Aril = -du A dR A d 2 uj . (7) 

where d 2 co = irii Am is the euclidian measure on the 2-sphere. 

Note. We have denoted by I, h, rh and rii the 1-forms l a dx a , h a dx a , rh a dx a and rh a dx a . We 
shall use this convention again. 

2.2 Rescaling of spin coefficients 

The rescaling of spin-coefficients under a general conformal rescaling 

~A t \ I A '•A A t \ 

g = \ l g , o = 1 1 o , i = i , oa = oa , la = * Ila , 

is described in [22] vol. 1 p. 359. The rescaled coefficients are obtained from the original ones 
by multiplication by a power of £1 with, for some coefficients, some additional terms that involve 
the derivatives of cu = log £1 along the original frame vectors. In the special case of the conformal 
rescaling (1), (6), £1 = R and these terms take the form 

Du = ! «V . = -i=(| + F|)(-logr) = -yfj ? , 
tu = ^V„„ = -i=(|-^)(-lo gr ) = 0, 



Sco = m a V a w = ^=(^ + -^-^-)(-logr) = 0. 

V2r sm ^ dtp ; 



D ' w - ' ! ° v " w = ^p(|- F l) ( - logr) = vf fi ' 
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and we have the following relations between the original and rescaled spin-coefficients 



k = r A K 


2 


7T = T7T 


p = r 2 p + \J\r 


a = ra 


A = A 


2 

a = r a 


(3 = rf3 


fi = p + \J~^R 


T = TT 


7 = 7"^ 


v = Rv 



The spin coefficients in the original tetrad (5) have been calculated in [19]. Using the array 
above, we obtain : 



A 



V = 7T 



A = o, 



m 



, 1 



5mR 2 - 2R 



cote 



-a 



2V2F ' ' 2V2F ' ' 2^/2 ' 

Note that the coefficients p and p were not zero for the original tetrad 3 . 



(8) 



2.3 Neighbourhood of spacelike infinity 

We work in the following domain for a given u$ « —1 

^uo := { ( n ' R > w ) ! u < u 0, < t < +oo , oj £ S 2 } . 

We foliate this neighbourhood of i° by the hypersurfaces (which are spacelike except for Hq 
which is null) 

% s = {u = —sr* ; u < uo} , < s < 1 . 

For s = 1, the hypersurface Hi is the part of the {t = 0} surface inside fi+ and for s = 0, Ho 
also denoted is the part of ,J? + inside ^„ . The level hypersurfaces of u within £l^ wm be 
denoted by S u , they are null. Given < si < S2 < 1, we will denote by Su 1 ' 82 the portion of S u 
between H S1 and H S2 . 

We need an identifying vector field between the hypersurfaces 7i s when decomposing 4- volume 
integrals over 0+ Q using the foliation. We use 

v = rlR 2 {l - 2mR)\u\- 1 d R . (9) 

It is tangent to the u = constant surfaces and is naturally associated to the parameter s in that 
u(s) = 1. The splitting of the 4-volume measure dVol 4 corresponding to the foliation {%s}o< s <i 
with identifying vector field v is the product of ds (being the measure along the integral lines of 
u a ) and i/jdVol 4 |^ s = r 2 R 2 (l — 2mi?)|u|~ 1 diid 2 w|% s (which is the resulting 3-volume measure 
on each % s ). 

We recall from [16] the controls we can infer on u, R and r* in the domain fi+ for |« | large 
enough. 

Lemma 2.1. Let e > 0, then for u < 0, |uo| large enough, in the domain Q+ o , we have 

r < r* < r(l + e) , 1 < Rr* < 1 + e , < <l + e, 1 - e < 1 - 2mi? < 1 , 

3 The property that p — was clear without calculation since p represents the geodesic expansion along the 
flow of I (which is a geodesic flow). This is clearly zero since for g the surface of the 2-spheres orthogonal to I and 
n is constant. As for fi, it is equal to — p', i.e. corresponds to the geodesic contraction along the flow of h. It is 
therefore also obviously zero for similar reasons. 
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and of course 



< s = — < 1. 



The factor r 2 R?{l — 2mR)\u\ 1 appearing in the expression of the vector field v satisfies 

<r 2 R 2 (l-2mR)\u\- 1 < . 



u u 



2.4 The Morawetz vector field 

The name "Morawetz vector field" is slightly inadequate. It refers to a vector field that is 
timelike in the neighbourhood of spacelike infinity and is transverse to ,y + . It is constructed 
from the actual Morawetz vector field in flat spacetime (see [17]) by expressing it in a coordinate 
system resembling our u, R, uj coordinates and brutally keeping the expression on the rescaled 
Schwarzschild spacetime. More precisely, the Morawetz vector field on Minkowski spacetime is 
defined by 

K = (r 2 + t 2 )d t + 2trd r 
and finds its simplest expression in the coordinates u = t — r,v = t + r: 

K = u 2 d u + v 2 d v . 

This is a conformal Killing vector field of Minkowski spacetime and is precisely Killing for the 
Minkowski metric n rescaled using the conformal factor Q = 1/r (i.e. f) = (l/r 2 )i]). If we use 
the coordinates u = t — r and R = 1/r, the vector K takes the form 

K = u 2 d u - 2(1 + uR)d R . 

We define the "Morawetz vector field" on the rescaled Schwarzschild spacetime in the coordinates 
u = t — r*, R = 1/r, as 

T a 8 a := u 2 8 u - 2(1 + uR)d R . (10) 
It has the following decomposition on the tetrad Z, n, m, rh : 

Ta = V? ( 2(1 + uR) + \ {uR)2p ) ia + u2 V? ha ■ 

Since K is Killing for the rescaled Minkowski metric R 2 n and since the Schwarzchild metric is 
asymptotically flat, the vector field T should provide an approximate Killing vector field, near i° 
and J^, for the rescaled metric g, which is precisely the reason why it was introduced in [16] to 
study the peeling of scalar fields on the Schwarzschild metric. A calculation of its Killing form 
shows that it is indeed a good approximation of a Killing vector in the vicinity of i° and J? : 

V (a T b) dx a dx b = 4mR 2 {3 + uR)du 2 = 8mi? 2 F _1 (3 + uR)lJ b dx a dx b . 

Note that a similar construction, based on the null coordinates (u = t — r* , v = t + r*) instead of 
(u,R), and also referred to as the Morawetz vector field, was used by Dafermos and Rodnianski 
in [7]. 
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3 Dirac and Maxwell fields 



A Dirac spinor field is the direct sum of a neutrino part \ A ' and an anti-neutrino part tpA ; in 
the massless case, the two parts decouple and Dirac's equation reduces to the Weyl anti-neutrino 
equation 

V aa 'tP a = 0. (12) 

Similarly, in the source-free case, the anti-self-dual part (f)AB = <P(AB) and the self-dual part 4>A'B' 
of the electromagnetic field decouple and Maxwell's equations are equivalent to the equations for 
4>AB ■ 

V AA > AB = 0. (13) 

Both equations are conformally invariant : spinor-valued distributions ipA and 4>ab = §{AB) 
satisfy respectively equations (12) and (13) on the exterior of the black hole if and only if the 
rescaled quantities ipA = Q~ 1 tpA = ripA and <$>ab = Q~ 1 4>AB = r 4'AB satisfy on the same domain 
the rescaled equations 

V AA 'i> A = 0, (14) 
V aa '4>ab = 0, (15) 

where V is the Levi-Civita connection for the rescaled metric g. 

3.1 Rescaling of the field components 

We decompose the physical Dirac field if) a and Maxwell field 4>ab onto the spin- frame {o A , l a } 
and the rescaled fields ipA = ripA and 4>ab = r 4 > AB onto the rescaled spin-frame {d A , t }. The 
decomposition is as follows : 

IpA = iplOA ~ ipOl^A , 

4>A = rip a = ripiOA ~ ripQLA 

= lpid A ~ tpO^A = i>\OA - TpoRiA , 

Hence, 

4>o = r 2 ip Q , V>i = ripi . (16) 

A simpler version is the following : 

ipo = tPao A = rip A ro A = r 2 ip , ^ = ^ A i A = ripA^ A = ripi . 
As for the Maxwell field : 

4>0 = r 3 4> , 4>i = r 2 (pi , (p 2 = r<j) 2 ■ (17) 

3.2 The rescaled equations 

The rescaled Weyl equation (14) can be expressed using the tetrad (3) and the associated spin- 
coefficients as follows (see S. Chandrasekhar [1]) : 

' b'^o - + (A- 7)^0 + (f - p)$i = o , 

< 

k Dtp! - 5'4> + (a - 7t)^ + (£ - P)ipi = , 
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the link being 

o = v AA 't A = (iy-k - Hi + (A - 7)^0 + (f - ~o A> 

+ (Dip! - S'tpo + (&- 7t)4j + (e - p)tp^j l A ' . 
This gives us the system ( with F = 1 — 2mR): 

% (d u + \R 2 Fd R ) ^--t=(d e + \ cot e + ^d v ) & - 

[ -yf^da^ -^ 2 {de + hote- ±- e d 9 ) & + ^1 = 0. 
This can be simplified as follows : 

b ip - 8-01 = , 

W>i- 8^0 = 0, 



(18) 



2V2F 



(19) 



where p, p , 9 and 5' are the weighted differential operators of the GHP formalism (Geroch-Held- 
Penrose [11], also referred to as compacted spin-coefficient formalism in Penrose and Rindler [22] 
Vol.1 section 4.12), which, applied to i/jq and tpi take the form 



R 2 F 



d R + 



2R - 5mR 2 



3^1 



V2 



1 i \ - - - 1 
+ - cot + — Vi , SVo = -7= 

2 sin 6* / ^/2 



d e + - cot 9 - -r—^dp Vo • 
2 smfe 1 



For the rescaled anti-self-dual Maxwell system, we have 

= V AA '<t>AB = (p'fo ~ Hi + (A - 27)<^o + 2f0i - <70 2 ) 5 A 'o B 

- (-D01 - <5' 0o + (2a - 7r)0o - 2p0i + £02^ o a 'l b 
+ - <50 2 - z>0o + 2A0i + (f - 2/3)0 2 ) Fo B 

- (l>0 2 - S'4>i + A0 O - 2tt0! + (2e - /5)0 2 ) F'tB , 
so equation (15) is equivalent to the system 

' D'4> - <50i + (A - 2 7 )0 O + 2f0i - o-02 = , 
£>0i - <5' 0o + (2d - tt)0o - 2p0i + k0 2 = , 
£>'0i - 502 - z>0o + 2A01 + (f - 2/3)02 = , 
k £)02 - 8' 0i + A0 O - 2^0i + (2e - A)02 = . 
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In the Geroch-Held-Penrose formalism, this takes on the simpler expression 





-901 


= 0, 


Ui 




= 0, 






= 0, 


\>4>2 


-g>i 


= o, 



(20) 



with 
t>'<£o 



-{d u + -R 2 Fd R )<fo -= — <p , d0i 



\>'ii = Jj(du + lR 2 Fd R )fa , H 2 



F ^ Tfl 



Or. 



1 



(de + -r—zd^fa , S'0o = -7= {do + cot — <9< 



V2 



sin 9 



V2' 



sin 9 



■><p)<P0 , 



d(f>2 = —i={de + cot 6* + -r—rd^fa , 3'</>i = —^=(0$ 



V2 



sm I 



sm ( 



3.3 Conserved quantity and estimates for perturbed equations 
3.3.1 The Weyl equation 

The conserved current for the Weyl equation is J a = tp A il) A , which gives the following closed 
3-form by contraction with the 4-volume measure (7) : 

u := *J a dx a = JjdVol 4 



td a + 4j h a d a - M m a d a - ip 1 i> rh a d a jdVol 



'4n 



= — \ipi\ 2 l A d 2 LO + IV'ol 2 "' A d 2 oo — i^iV^o I A h A m + i^oV'i I A n A m 



di? A d 2 w 



R 2 



^0 



dti A d 2 w 



+a/23? (vWi) dn A di? A sin 0d</? + \/23? (iipo^j du A di? A d9 . 
On a given hypersurface ~H S , we have 



(21) 



di? = — du 



and the conserved quantity takes the form 



- -1 ) i? z |^ | + 



^ - 1 ] lr 

\u\ 





2 






V>0 


+ 




V 



— du d 2 w 
2 

' "^"dw d 2 w . 



(22) 
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On S u , 



'Su 



di?d 2 w. 



(23) 



Lemma 3.1. The energies on S u , u < uq and % s , < s < 1 have the following simpler 
equivalents (meaning that there are constants independent of u < uq, < s < 1 and the smooth 
spinor field ipA such that the energies on S u and T{ s are controlled above and below by these 
constants times the simpler expressions) : 



dRd 2 u, 



R 



du d 2 u! . 



(24) 
(25) 



l]-oo,u [ u xSZ \\ u \ 

Proof. This is a direct consequence of lemma 2.1 and of the fact that 

1 2 2 
_<__!<_. □ 

s s s 

The closedness of the 3- form u> gives for any smooth solution ipA of (14) with compactly 
supported initial data : 

£n n (VO + £ s %>2 (j>) = £n S2 (YO for any < si < s 2 < 1 . (26) 

Now, consider a Dirac equation with error terms of two types, a potential P and a source Q : 

V AA 'i, A = P AA '^ A + Q A ' . (27) 

Then, differentiating the current 1-form J a , we get 

V AA ' [Ma) = (p aa> + P AA ) + 2R (q a '^a) 

which is an approximate conservation law. When integrating this over the 4-volume flti ,S2 
bounded by hypersurfaces H Sl , H S2 , f° r < si < S2 < 1, and the part Sul' S2 of S UQ between % 
and / H S , 2 , we obtain for any smooth solution of (27) with compactly supported initial data : 



Sl 



(( 



pAA> + pAA^J A> + m (^qA'I 



dVol 4 



< (1 + e) 2 jH | (P AA ' + P AA ') Ma> + (q a '^ a ) 



-dud uids . 



Energy estimates will be established using the Gronwall inequality provided the integrand on 
the right-hand side can be estimated by the energy density on % s , in a sufficiently uniform way 
so as not to prevent integrability in s. 

3.3.2 The Maxwell system 

An anti-self-dual Maxwell field 4>ab has a stress-energy tensor given by the following expression 

Tab = 4>ab4>A'B> ■ 
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In order to define an energy current, we need to choose a timelike vector field to contract the 
stress-energy tensor with. A natural timelike vector field would be the Killing vector d u which is 
equal to dt in the Schwarzschild coordinate system, but d u becomes null on and we require 
more control there in order to establish peeling results. So we use the Morawetz vector field (10) 

T a 8 a = u 2 8 u - 2(1 + uR)d R . 

The associated energy current is the vector field V a = T o6 T& whose decomposition on the 
Newman-Penrose tetrad l,h,rh,rh is given by 

V a = ^ll(2 + 2uR + ^F)\^ + ^lu^ 2 \ 



- ( K l-{2 + 2uR + ^^F)\M Z + \j^\^jh a 
| (2 + 2uR + ^^F)^J + ^ u %J^j m a 

j{2 + 2uR + ^^F)4> J 1 + ^tffa'faj m a . 

Since T a is not an exact Killing vector, V a is not divergence free and it satisfies merely an 
approximate conservation law 

V a V a = V {a T b) T ab = 8mR 2 F~ 1 (3 + uR)T ab l a l b = 8mR 2 F~ 1 {? > + uR)\fa\ 2 . (28) 
The energy 3-form is the Hodge dual of the energy current 
u := *(V a dx a ) = VjdVol 4 

(2 + 2uR + ^-F)\k? + ^1 u 2 \H^ HAdM 

+ [ \/ ^(2 + 2uR + { -^F)\M 2 + yjl u 2 ^ 2 ) (h A d 2 ^) 



f ! \/ -(2 + 2uR + -^-^—F)(j)i(/) + y -u'fafaj (-U An Am) 

+ (^fji 2 + 2u R + ^Y-F)4>o4>i + y|v0i0 2 ^ (if Ah Am) , 

(recall that d 2 w = im A m is the euclidian measure on S 2 ). 
On a u =constant hypersurface S u , the energy is given by 

SsM--= jf w = jf ^(2 + 2^ + ^F)| < Ao| 2 + y|n 2 | ( A 1 |^di?d^. 
On a u = — sr* hypersurface % s , recall that 

di? = —clu 



(29) 
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and therefore (for < s < 1) 



+ [2 + 2uR+ ( uR ) 2F ) |^| 2 + -n 2 |0 2 nd U d^. (30) 



As s — )■ 0, this expression simplifies to give the energy on 

S J+ Q>):={ w=/ ^f + ^l^dudV (31) 



u u 

using the fact that 

(uR) 2 



s 

?2 , 



= {—uR)Rr* — > as r — >■ +oo with u bounded 



and similarly R 2 /s — > as r — > +oo. 

Lemma 3.2. We have the following equivalent simpler expressions fo the energies on S u , u < uq 
and % s , < s < 1 

£sM - [ (|^o| 2 + ^i| 2 )diM 2 ^, (32) 
SnM ^ (^|0o| 2 + |0i| 2 + « 2 |02| 2 )d U d 2 a;. (33) 

Proof. We notice that 2 + 2uR + Mlf-F vanishes for 

2 , 

uR = - — (l±V2mR). 

r 

We know from lemma 2.1 that in ^„ , — 1 — e < uR < with < e « 1 and also 



2(1 - e) < -(1 ± V2mR) < 2(1 + e) . 
r 



It follows that 2 + 2uR + v 2 ' F vanishes nowhere for < s < 1 and tends to 2 as uR — > 0, so 
this quantity is also bounded below away from zero and above, uniformly on f2+ Q . The lemma 
then follows from lemma 2.1. □ 



4 Peeling 

4.1 Peeling for Dirac 

We have already established energy estimates for ipA between J?^, S Uo and T-L\. Now we establish 
estimates for successive derivatives of ipA ■ We do not need to commute all directional derivatives 
into the equation ; as our goal is to control transverse regularity on ,J? + we focus on derivatives 
in the direction of Or, i.e. of l a . We denote by Dr^a the spinor 

Dri/ja ■■= dRipiOA - OriPola 
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and we work out the equation satisfied by D R ip. In order to obtain a more useable expression 
for this equation we multiply the first line of (19) by y/2F, keep the second as it is and commute 
d R into the resulting system. We get : 



( (2d u + R 2 Fd R ) d R ip - y/wtdnfa + R(1 



5mR 



) d R 4i 



%(d R - ft) d R ^ - &djnh 

This can be re-written as 



= -2R (1 - 3mR) d R ^ - m^|.3^i - (1 - 5mR) 4> 

W - ( 2 ^)3/2 ^1 • 



/2F 



2R~hmR 2 



9 R TpQ - dd R 4>i 



4,R(1 - 3mR)d R ^ ~ Wi 



(d R - d R ^ - W8 R ^ 
and as a spinorial equation takes the form 



l-5m.fi 
V2F 



1>1 



V'O, 



V AA '(D R ^ A ) = -(^R(l-3rnR)d R ^ + ^Ui + ^^^]o 



m 



d R tjJi + 



m 



rV'i t 



/2F (2F)3/ 2 
The conservation law associated with equation (34) is the following 

V AA ' [(d r ^ [d r ^ a ^\ = 2K [(v aa 'D r ^a) D r ^ a 



(34) 



2K 
+ 



2 /v 777 

-R{1 - 3mR)d R ^ + j^i - 



1 — 5mR 



2F 



V'O d R if) 



m 



d R 4>i + 



m 



(2F) 3 / 2 



ipi d R ipi 



In order to obtain energy estimates for D R ip using Gronwall's inequality, we need to estimate the 
right-hand side by the energy densities for ip or D R ip. Two types of terms present a difficulty : 
those involving angular derivatives and the term involving ipQ and d R ^Q without a factor of R. 
In order to control the angular terms, we must commute angular derivatives into the equation ; 
this will turn out to give us the additional control on ip that we need. 

The commutation of angular derivatives into the Weyl equation is best described, and per- 
formed, using the GHP formalism, i.e. we use the form (19) of the Weyl equation. In our 
framework, using the values of the spin coefficients, we have the following identities 

p , 8 7] = for rj of weight { — 1, 0} and p , 6' rj = for rj of weight {1, 0} . (35) 

So when we apply 9' to the first equation and 9 to the second, we obtain 



ps^i - m'tpo = o . 



(36) 



14 



Note that a^o has weight {0, 1} and 6^1 weight {0,-1}, so putting 

(D u iP)a> = S^lOA' - , 
equation (36) can be written as a conjugate Weyl equation 

V AA '(Dj) A ,=0, 

i.e. is a symmetry operator for the Dirac equation on the rescaled Schwarzchild metric, 
sending anti-neutrino fields to neutrino fields. 

Now consider some smooth solution tpA of (14) with compactly supported initial data. Both 
V> and D^ip satisfy the energy equality (26), which entails for any s £ [0, 1] : 



j n (j^l 2 + l^i| 2 )d«d 2 w ~ jf ^|^,| 2 + |^i| 2 )d«d : 



U) 



+ 



hi (l«l 



Wo 



dRd 2 oj , 



R 



8'^ | 2 + |6^i| 2 ) dnd 2 w 



R 



Us \\ u 



8'V> | 2 + 1 1 2 ) dud 2 w 



+ 



3^o 



di?d 2 w . 



Note that this immediately gives us a control of the 4-volume L 2 norm of i/jo and S'^o m terms 
of the energy on "Hi using the foliation by S u with the identifying vector field d u . However, 
this works only once we have obtained the estimates, it does not allow to control a perturbed 
equation with an error term (^o) 2 i n the conservation law. 

Integrating on Qu ' S2 the conservation law associated with equation (34), we obtain 



< 2 



S2 

si Jn, 
+ 



%R{1 - 3mR)d R i> - ^8^i - 1 ^ R 4> 



m 



'IF 



m 



2FV2F 



\9r^i\ ) i — r dud 2 u;ds 

' \u\ 



The last two error terms are trivially controlled by the energies of ipA an d Dri^a j the first error 
term, thanks to the l/\u\ coming from the Leray measure, is exactly controlled by the energy of 
Dri^a- The difficulties are with the second and third error terms. For the second term, we use 
the control we have obtained over angular derivatives as follows 



si Jn 



\dRipo\-^-dud 2 ujds < 
\u\ 



< 



< 



r 2 i r 
J si V$ Jn 

J 8! Vs Jn 

i v 7 ^ 



\dRipo\\ -. — rdnd cods 
u\ 



R 



¥i| + t-tI^oI dud 2 ujds 
\u\ 



2 -1= ( £h, (A, Y>) + £ Ha (Dr4>) ) ds , 
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which allows to apply a Gronwall inequality since 1/y/s is integrable on [0, 1]. The third term is 
the trickiest. We use the fact that the lowest eigenvalue of S' on weighted scalar fields of weight 
{1,0} is positive (see [22] section 4.15) ; this implies 



/ |^ | 2 d 2 ^< / |6'V> | 2 d 2 ^ 
Js 2 Js 2 



is 2 Js 2 

uniformly on 0+ whence using the same method as for the previous term, 



«(}' 

r s 2 



/ / \d R ^ \--dud 2 u;ds < / -= / (SVo + ^{dR^dud^ds 

J S1 Jus V2F \u\ J ai J Hs \u\ 

and since 9'^0 = t^i which is controlled uniformly in fi+ by |<9ijY>i| + |V>i|, we can once again 
apply Gronwall's inequality. Successive applications of d R will produce error terms which can 
be controlled by similar techniques using angular derivatives and lower order norms. We get the 
following result : 

Theorem 1. There exist positive constants C n , n £ N such that for any smooth compactly 
supported data on ~Hi, the associated solution t/m of Dime's equation satisfies 

n p n p 



E E £ K, ( D tD P R %A) < a E E «Ki ^ D t D R k ^A) , 

p=0 k=0 p=0 k=0 

n p 

£ £ ^ (jWSko < c?n E E ( D Z D R k M + ^ (^r^)) 



p=0 fc=0 p=0 fc=0 

T/iis extends to the spaces of initial data rj n ("Hi) obtained by completion o/Cq°(Hi) wi norms 

\ 1/2 



p=0 k=0 



4.2 Peeling for Maxwell 



Contrary to the case of the Dirac equation, we already have a problem with the basic energy 
estimate since the energy current V does not satisfy an exact conservation law. However, the 
error term in the approximate conservation law satisfied by V is easily controlled by the energy 
density on the hyp ersur faces 7i s . This gives us the following result : 

Proposition 4.1. There exists a positive positive constant C , such that for any smooth data 
on Hi satisfying the constraints and supported away from i°, the associated solution 4>ab of 
Maxwell's equations satisfies for any < s < 1 



£h s {<Pab) < C£ n i((pAB) , 
£-Hi{4>Ab) < C (£-Hs{4>AB) + £ s s,i((f>AB)^ 



We now obtain similar estimates for successive derivatives of the Maxwell field, starting with 
the derivative with respect to R. We multiply by y/2F the first and third equations of (20) and 
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commute d R into the system. We obtain 



f {2d u + R 2 Fd R - (5mR 2 - 2R)) d R fo - y/2FZd R fa 

= -2R(1 - 3mR)d R o ~ 2(1 - 5mR)4> 







'R9l 



(2d u + R 2 Fd R )d R (pi - y/2F&d R fa = -2R(1 - 3mR)d R fa - J |mcV 2 , 



{ -y/$d 2 R <h + jlpd R <h-&d R <f> 1 
which we rewrite as 

\>d R (i)o - 6<9i#i = 



d 



RH>2 



(2F)3/ : 



■4>2 , 



j,R(l - 3mR)d R( f>o - Jf(l - 5mR)cf> - fdfa , 



Vdnfa - hd R 4> 2 = -V2FR(1 - 3mR)d R 4> l - fBfc , 



Putting 



Dr^AB ■= d R (f)Qi A i B ~ d R <pi{d A iB + tAOB) + d R (p 2 OAOB , 
the above system is the perturbed Maxwell equation for D R <pAB '■ 

V AA 'D R j>AB = -(y2FR(l-3mR)d R fa + jdfo)d A 'd B 

+ (JjRQ ~ ZrnR)d R ^ + - 5mR)0 o + ^) 6 A 't B 



The associated approximate conservation law is therefore 

V A4' (t bb '(D r ab )(dJ a , b ,)) = V {a T b \D R j> AB )(Dj A , B ,) 



+2K {T BB ' [- (V2FR{\ - 3mR)d R 4> 1 + ^3<£ 2 ) 6 A ' 



OB 



+ f J jR{l - 3mR)d R 4> + A/|(l - 5mR)4>o + jUi ) d A 't B 



m 



2F 



(2Ff/ 2 



(j/.x > 2 f — — 7-702 1 i OB + —=(>l!0:l L B 



D R 4> A > 



17 



which in full details, using the decomposition (11) of T a on the rescaled tetrad, reads 



V AA> (t bb '(D r ab )(D r } a , b ,)) =8mR 2 F~ L (3 + uR)\d R <j)o\ 



~ ^2(1 + uR) + ( uK ? F ^ (1 _ 3mR)R\d R 4> \ 2 

- (2F 3 / 2 (1 - 3mR)\Ru\\u\ + j (2(1 + uR) + 

-mu 2 \d R (j)2\ 2 
-23? 



|^i| 2 



u 2 d4> 2 d R 4>i + ^u 2 4> 2 d R 4>2 



+ 1(1 _ 5 mi?) ^2(1 + uR) + M>!^ fo d J Q 

V2m ( , „, (uR) 2 F\ ~ 7 n ~ 
+^ (2(1 + + J S^i^o 



(37) 



As was the case for the Weyl equation, it is necessary to control the angular derivatives first 
before gaining a control on the derivative with respect to R. However, the weights of (fro, (\>\ and 
4>2 are {2,0}, {0,0} and {—2,0} and for an anti-self-dual field <\>a>b' we would get the weights 
{0, 2}, {0, 0} and {0, —2} ; since the weights of 6 and 6' are {1, —1} and { — 1, 1}, by commuting 
these operators into the Maxwell system we do not recover any of the adequate weights of the 
components. So we cannot hope to commute a well-chosen arrangement of 9 and 3' into the (anti- 
self-dual) Maxwell system and obtain another (self-dual) Maxwell system. But the Schwarzschild 
metric is spherically symmetric, we have a 3-dimensional space of Killing vector fields tangent to 
the sphere which generate all rotations. We choose a basis of this space, denoted X, Y, Z, such 
that —X 2 — Y 2 — Z 2 is controlled below and above by the positive Laplacian on S 2 . Since X, Y 
and Z will commute with Maxwell's equations, this gives us a control analogous to proposition 
4.1 over angular derivatives of any order. 

Integrating the error term of (37) over ^ {r < s < 1}, for < r < 1, and splitting 
the result as an integral in s of integrals over the hypersurfaces ~H S , we gain a factor 1/u in the 
3-volume measure. With this, the integrals of all terms on % s are controlled for each s by the 
sum 

£h s {4>ab) + £h 3 {D R (Pab) + £h s {Dx4>ab) + £u s {D y <Pab) + £n s {Dz4>ab) , 

except for the last two terms. The treatment for both terms is similar to the Dirac case. For the 
first, we use the fact that the L 2 norm of c/)q on the 2-sphere is (uniformly in f^„ ) controlled 4 by 
that of S'^o and then the second equation of the Maxwell system giving us the equality of (5'<fio 



4 The first eigenvalue of a on weighted scalars of weight {2,0} is again positive, see [22] section 4.15. 
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and So we obtain for < s < 1 



|(i _ 5mi? ) ^2(1 + uR) + M>!^ fo dR l 
1 f ~ R 

V* J-Hs V M 



- — -dnd 2 w 



< 



u.. 



iSVofdudW / |5^o| 2 7^dud 2 w 

\u\ 



< 



HI 2 + |9^i| 2 ) dudW J \d R ^\ 2 ^dud 2 oj 
< (Sh 3 (4>ab) +£h s {Dr4>ab)) ■ 



The second is simpler to deal with : 

\/2m 



2(l+«i2) + 



( Uj R) 2 F 



< 



Bcf)idR(f)o 
/ i? 

'l^ij^oU / 1 — :dud 2 u> 



i — rdud 2 u> 

\u\ 



4= / |30i 

~ ( d x4>ab) + £n s (Dy4>ab) + £w> (Dz<Pab) + £u s (Dr^ab)^ ■ 

Remark 4.1. ./Vote that since <f>\ is the component of <f>AB whose weight is {0,0}, <5(j)\ is merely 
5(f>i = Vtfj^ij it involves no spin- coefficient. And since m is a bounded vector field on S 2 , the 
L 2 norm officii on S 2 is controlled (uniformly on Q^o) by the sum of the L 2 norms of Dx4>i, 
Dyfpi and Dz4>\- 

The successive derivatives with respect to R will be controlled in a similar way, controlling 
first the angular derivatives. This gives the theorem : 

Theorem 2. There exist positive constants C n , n £ N such that for any smooth data on TL\ 
compactly supported, the associated solution <f>AB of Maxwell's equations satisfies 

£ Sj+ (D%D*D*D%$ab) <C n Yl ^ (D^D^D^D^ab) , 

k\+k2+kz+ki<n k\+k2+k^+k4 L <n 

X; £h, (D^D^D^D^ab) <c n Y, { £ K (d^d^dIw^ab) 

ki+k2+ks+k4<n ki+k2+k3+kn<n 

+£s U0 (D k x WpD^D k ^ 

This extends to the spaces of initial data f) n ("Hi) obtained by completion ofC^iTLi) in the norms 

/ \ V2 

W&abW^chi) = 



\ki+k2+k3+k4,<n 
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5 Interpretation 



In this section, we check that for a given order of transverse regularity at J^ + , our classes of 
data ensuring that the rescaled solution has at least this regularity are not smaller than they 
are in the flat case when the full embedding in the Einstein cylinder is used. Of course, since 
we have only worked in a neighbourhood of i° in Schwarzschild, what we really mean by this is 
a comparison of the asymptotic constraints on the fall-off of initial data. In the case of Dirac 
and Maxwell fields, this comparison is made easier than for the wave equation (see [16] for 
the analogous interpretation for the wave equation) because the energy on a spacelike slice is 
conformally invariant. 

Let us consider a 4-dimensional globally hyperbolic (and therefore admitting a spin-structure) 
spacetime (M,g), S a spacelike hypersurface, u a its future-oriented unit normal vector field, if)A 
a solution to (12) and <&ab a solution to (13). Let 0, be a smooth positive function on M, and 
put 

g := £l 2 g , ipA ■= ^~Va , 4>AB ■= ^ 1 4>ab ■ 
The unit normal to £ for g is now 

v a = fT V 

and if we denote by fi (resp. fx) the measure induced on S by g (resp. g), then 

fi = Q?ijl . 

The energy of the rescaled Weyl field on £ is given by (in this section, we shall denote by £ the 
energies for the rescaled metric and £ the energies for the unrescaled metric) 

£e(Y>) = / ^ipAtpA'dfi = / u a tp A tpA'dn = £ S (V0 ■ 
Jt, Jt, 

For the Maxwell field, we need a choice of observer (or merely of timelike vector field) t a in the 
neighbourhood of £ to define the energy and then 

= / l> a t b 4>AB<i>A>B>dfl= / V a t b (j)AB4>A'B'^ = £v(<P) ■ 

Jt, Jt, 
Note that the observer is not rescaled. 



5.1 Dirac fields 

We now compare the classes of data for Dirac fields in the flat case and in the Schwarzschild 
spacetime. In the flat case (m = 0), the conformal embedding of Minkowski spacetime into the 
Einstein cylinder is realized using the conformal factor 

9 

Q 



Vl + (t + r)Vl+(t-r) 2 

and the vector field used for increasing the regularity in the energy estimates is the time trans- 
lation along the Einstein cylinder 

d 1 / ' o 9 . d d 
87 = 2 + < +r) B t +2tr 3-r 
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This is a Killing vector field on the Einstein cylinder, i.e. a conformal Killing vector field of 
Minkowski spacetime. Let ipA be a Weyl field on Minkowski spacetime, denoting by £ the t = 
hypersurface and choosing the spin-frame given by the choice of Newman-Penrose tetrad 

1 1 1 i 

I = -^=(<9 t + dr), n = -j=(d t + d r ) , m = ^=(#0 + ]^ d v) > ( 38 ) 

the future-oriented unit (for the unrescaled metric) normal to £ is dt = -^(l + n) and we have 

= £sW0 = iy a ^A'd 3 x = ±= jj\M 2 + |^i| 2 )d 3 x . 
On the Einstein cylinder, we have the following energy equality 

and commuting d\ into the equation we get 

£x(d k J) = £Md k J) ■ 

The right-hand side is a measure of transverse regularity at ,J? + . 

We work out the constraint on initial data corresponding to the first level of regularity. 
Using the expression of d T in terms of (t, r) variables at t = 0, the left-hand side for k = 1 can 
be rewritten as 

1 + r 2 - 1 + r 2 

Using the ellipticity of the spacelike part of the Dirac equation (i.e. using the Bochner-Lichnerowicz- 
Weitzenbock formula on R 3 ), this corresponds, modulo lower order terms, to an L 2 control over 

(1 + r 2 )d r V> and ^^-V S 2ip 
r 

independently. This is what defines the scale of weighted Sobolev spaces obtained on S by 
requiring the finiteness of the energies of d^'ij). 

Note that this is not the original conformal peeling construction by Penrose, but it is the 
closest equivalent in terms of Sobolev spaces (see [16] for a more detailed presentation for the 
wave equation). 

In the Schwarzschild case, the energy equality is now 

and to raise the regularity up to order k, we use all combinations up to order k of 8r and D^. 
At first order, the quantity defining the space of data on "H 1 is 

£ n i(4>) + £ n i(d R 4>)+£ n i(Dj), (39) 

where D u can be replaced by V52, the former acting on the components, the latter on the full 
spinor. First of all, using the conformal invariance of the energy, for the unrescaled spin-frame 
associated with the tetrad (5), 

£ n i(4>)= [ (IVoI 2 + |^i| 2 )i 7ll/2 r 2 sin6'drd6'd( / 9 , 

Jn 1 
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which gives a control equivalent to the flat I? norm. Now up to lower order terms, using the 
conformal invariance of the energy, (39) is equivalent to the same expression with unrescaled 
energies for the unrescaled field ip. The last term in (39) bounds the I? norm of V s 2 ^ with 
no weight, which is weaker than in the flat case, but the second term will impose a constraint 
similar to the flat case, though in a more mixed manner. Indeed, using the equation satisfied by 
tp in components (see for example [19]) 

/ _i o \ 1 F' F 1 / 2 



we find that 



8r%1) = ^(d t + drjifj = ^ ® r *^ J ° J _ r p l l 2 D< JJ ip + lower order terms, 

J ° J ~ rF^ x l 2 D^ + lower order terms. 

We see that in terms of weights, this is similar to the flat case. The gain is that we only control 
r 2 d r applied to ipo and not to ipi. This is not surprising since tpo is the component of ij) that 
propagates dominantly to the left (towards the black hole) and which, were its fall-off at infinity 
not strong enough, might therefore propagate singularities along ,y + from i . Also the weight 
in front of the angular derivatives in the expression of 8r^ is the same as in the flat case, but 
the term is not controlled on its own, only in combination with r 2 d r ipQ. This is what defines the 
scale of weighted Sobolev spaces obtained from our energy estimates in the Schwarzschild case. 

The conclusion is that the constraints for peeling in the Schwarzschild case are weaker than 
the ones obtained in Minkowski using the full conformal embedding in the Einstein cylinder. In 
other words, peeling in Schwarzschild at any order is valid for a class of data slightly larger than 
the usual class in the flat case. 

Remark 5.1. This does not mean that the asymptotic structure of Schwarzschild allows peeling 
for more general data than in Minkowski. This only shows that our definition of peeling is more 
general than the one usually considered. The uniformity of our norms in the mass in any compact 
interval [0, M\ shows that for our definition, the classes in Minkowski and Schwarzschild are the 
same. 



5.2 Maxwell fields 

In the Maxwell case, the essential ingredients are the same as in the Dirac case : conformal 
invariance of the norm, same vector fields used to raise the regularity and controlling better the 
part of the data propagating to the left, plus one more ingredient, the Morawetz vector field. 
What is not clear for Maxwell fields is that the basic energies should give equivalent controls. 
The reason why this is true is that in flat spacetime, the Morawetz vector field and the time 
translation along the Einstein cylinder differ by a constant multiple of dt (see [16]) which gives 
a weaker norm than the two others, plus the local equivalence of the norms in the mass m. Let 
us describe this more explicitely. 

On the Einstein cylinder, we use the vector field d T both for defining the conserved current 
for the Maxwell field and for raising the regularity in the energy equalities. Using the conformal 
invariance of the energy and the expression of d T , decomposing the unrescaled field in the spin- 
frame associated with (38), we get 

= / (|0 o | 2 + 2|0i| 2 + |0 2 | 2 )^^d 3 x. 
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In the Schwarzschild case, we define the energy using the Morawetz vector field given in terms 
of the vectors of the rescaled tetrad by (11). The initial energy therefore is given by (using the 
fact that the future unit normal to the hypersurface {t = 0} for the metric g is -^=(1 + n) and 

also u = — r* on H 1 ) 

£ H ^) = £ H ,{<t>) = ^ i (^(2(l-r^) + ^|^)r 2 |0 o | 2 



which in S7+ Q is equivalent to 

I (|<Ao| 2 + |0i| 2 + |<A 2 | 2 )r 4 drdw 
Jn 1 

and therefore gives the same control as the flat norm involving d T instead of the Morawetz 
vector field. When raising the regularity in the energy estimates, the principles are the same as 
for Dirac. On Minkowski spacetime, a dt applied to (j) gives a control on the full gradient on R 3 
because Maxwell equations entail that the d'Alembertian of the field vanishes and this ensures 
the ellipticity of the spacelike part. So the flat spacetime energy at the first level of regularity 
satisfies 

£*{<!>) + Svi^-dtt) ~ J (H 2 + (1 + r 2 )|V R3 0| 2 ) (1 + r 2 )d 3 * . 
In the Schwarzschild case, the first level of regularity is controlled by the energy 

£ n i(j>) + i nl (d R j>) + s n i(Vs4), 

which is equivalent to 

Similarly to the Dirac case, we can evaluate d^cp using 

Or = —(dt + dr.) 

and Maxwell's equations. We have : 

9r4>o = — (V 2Fm a d a (f>i + 2d Tf (f)o) + lower order terms, 
F 

9r4>i = d-r m a d a (j)Q + lower order terms, 

9r4>2 = \j^r 2 m a d a 4>i + lower order terms. 

So we observe exactly the same phenomenon as for Dirac, a similar control as in the flat case, 
except that the control on radial derivatives is weaker : only the component propagating to the 
left is explicitely controlled on the initial hypersurface. 
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5.3 A remark on the constraints 

This paper is concerned about the asymptotic assumptions on initial data ensuring peeling at a 
certain order and on comparing these asumptions in Schwarzschild and in flat spacetime. There 
is one aspect, for Maxwell fields, which may affect the classes of data leading to a peeling at 
a certain order in different ways in Schwarzschild and Minkowski : the constraints. We do not 
address this question in detail here but we simply make a remark on the integrability of the 
constraints from infinity in a reasonably large class of functions and the compatibility between 
the constraints and our energy norms. Let us consider the Maxwell system in the unrescaled 
Schwarzschild spacetime using the tetrad (5). The energy at t = associated with the Morawetz 
vector field in the r* > — no region is equivalent to 

/ (|0o| 2 + |0i| 2 + |<A2| 2 )r 4 dr,d 2 w. 

Taking the difference of the second and third equations of the system, we get 




rescaling the field components by r 2 , 

fa := r 2 4>i , i = 0,1,2, 

this reads 

Fd r fa = - t (902 - g'0 O ) • (40) 

Decomposing the field into spin-weighted spherical harmonics, the operators 3 and 8' will turn 
into multiplication operators by constant factors. If we consider a simplified situation where 
the behaviour of the field at infinity is in powers of r _1 , the finiteness of the Morawetz energy 
implies that the rescaled field components fall off at least like r _1 . Assuming such a behaviour 
for <po and an d integrating (40) from infinity starting from the value zero, we recover the same 
fall-off for 4>\. 

So heuristically, the constraints are compatible with the Morawetz energy in both the Schwarz- 
schild and Minkowski spacetimes. This heuristic argument does not address the question of poly- 
homogeneous solutions of the constraints. Our sole purpose here is to show that in Schwarzschild 
and Minkowski, there will be large classes of solutions of the constraints compatible with our 
function spaces and that the constraints will not introduce in the Schwarzschild case any addi- 
tional restriction compared to the flat case. 

A Covariant derivative approach 
A.l Curvature spinors 

We calculate the curvature spinors for the rescaled metric. Recall that given a spacetime (M,g) 
with a spin structure and equipped with the Levi-Civitta connection, the Riemann tensor R a bcd 
can be decomposed as follows (see [22]) : 

Rabcd = XabCD £A'B'£C'D' + &ABCD' £A>B'^CD + ^A'B'CD ^AB^C'D' + Xa'B'C'D' ^AB^CD , 
where Xabcd is a complete contraction of the Riemann tensor in its primed spinor indices 

Xabcd = -^Rabcd£ A B £ C D 
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and $ a 6 = <J>( a 5) is the trace- free part of the Ricci tensor multiplied by —1/2 : 

2$ afe = 6Ag ah -Rat, A = ^Scal 9 . 
It is convenient when dealing with conformal rescalings to use instead of <3? a b the curvature 2- form 

Pab = $ab ~ kg ab 

because of its simpler transformation law. Also, it is usual to isolate the totally symmetric part 
^ABCD of Xabcd, referred to as the Weyl spinor, which describes the conformally invariant 
part of the curvature : 

Xabcd = ^abcd + A [eac^bd + £ad^bc) , ^abcd = X^ A bcd) ■ 

The scalar A and the curvature spinors P and ^ have simple rules of transformation under a 
conformal rescaling g = 0, 2 g, given by (see [22] p. 120-123) : 

^ABCD = ^ABCD , 

A = fr 2 A + jQ~ 3 an , □ = v a v a , 
P ab = Pab - v 6 T a + r AB ,r BA > , with r a = n^Vofi = v a log n . 

Lemma A.l. For the rescaled Schwarzschild metric (1), the values of A and are : 

A = mR/2, 
$ ab dx a dx b = 1 ~^ mR (R 2 Fdu 2 - 2dudFc + du 2 ) 

= (1 - 3mR) (hg + da; 2 

Proof. The value of A was calculated in [16]. In order to evaluate 4> a fe, we determine P a b. 
First note that the Schwarzschild metric is Ricci flat, whence = P a f, = R a f, = 0. Hence, 

P a b = -V fe T a + Yab'^ba' ■ 

Since Q = R = 1/r, 

dr 

T a dx a = --. 

r 

We need to determine its spinor components. We do so in the dyad {o A ,i A }, denoting x° = t, 

x 1 = r, x 2 = 9, x 3 = ip : 

1 11 BB 1 ^ BB' 

lAA' = 9AA> = 9 £AB£A'B>9i = — £ AB£A'B'9i 



and 



It follows that 



BB' = ( n i ~ m i \ 1 / I 
91 ~\- mi h )~ V2F V - 1 

1 [F ( 1 



Y AA 



-1 
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The non zero components of a ab := Yab'^BA' are therefore 

F 

"OO'OO' = ^00'"^00' = Oll'll' = Tn/Tii/ = > 

F 

"oi'io' = Too'Tn/ = aio'oi' = Tii/ Too' = — ■ 
and the 2-form a a & reads : 

Tab'Tba'cIx'Mx 6 = — (I a l b + n a n b - m a fh b -m a m b ) dx a dx b 



F 
~2 

i? 2 F 2 , o „ , , „ dR 2 F 



R 2 lJb + R 2 n a fi b — rh a m b — m a rhbj dx a dx b 



du 2 - FdudR + — 5= do; 2 . 



2 R 2 2 



We now calculate VbT a : 



V,,T ( ,d./-"d - \:, I dx 6 = ^ + iri b dx a dx b 



and among the Christoffel symbols 

r i _ 1 ic ( ^ffac dg hc dg ah \ _ F / d gal d gb i 9g ab 
ab 2 9 \dx h dx a dx c )~ 2 \8x h 8x a Or 



the non-zero ones are 



Hence 



FF 1 f 1 
r 00 = — > = -— , r^ 2 = -Fr, r* 3 = -Frsin 2 . 



/? 7?' /?' 

V 6 T a dx a dx fe = -=- + dt 2 dr 2 -Fdw 2 

r 2 2r 2rF 



^ + 5 ~ F- 1 ^ 2 ) ~ Fdu 2 

dR 2 

FmR^du 2 - 2mRdudR + —5- - Fdu 2 . 

R z 



From this we can infer the value of & ab : 

$> ab dx a dx b = P ab dx a dx b + Ag 

„ „ „ mi? „ 
= T^b/Tb^/ - VfeTa H —5 

= 1 -* mR (R 2 Fdu 2 - 2dudR + du, 2 ) . 

This concludes the proof. □ 

A. 2 Energy estimate for the transverse derivative of a Dirac field 

First we set up a general formula for the commutation of a directional covariant derivative into 
the Dirac equation, then we apply it to the derivative along 8r. 
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Consider any vector field V a and commute Vy into the Weyl equation : 

= V a V a V BB '^ B = -V a V a i B ' c 'v C c4 C 

= -i B ' C ' (v a A a J c + Vcc> {y a V a ^ c ) - (il c V a ) V a ^ c ) 

= -i B ' c 'v a A ac i, c + v BB ' (v v 4, B ) - (vVj Vj B 

where A a ^ = V a V& — V 6 V a . The first and third terms in the right-hand side can be calculated 
more explicitely : 

V a A a J c = V a 



£A'C'X ACE +EcA®a>C'E 



The symmetries of the Riemann tensor imply that 



v c OA- 3mR „ 
■X-ace = SAeae = — 2 — £j4E ' 

whence 

-e V ea'c^ace V ~ £ —^—Vc'VA-^—V ip A - 
The term involving can be written 

-e B,C "iAC*A>C>E C 'ip E = -VcA^ EB ' CA '^E = -VM B '^B • 
It follows that the equation satisfied by Vyip B is 

v BB ' (v v iPb) = (vVj + v a ¥ b i> B - ^v b i> B . (4i) 

If the vector field is V a d a = d R , then 

V a y b = g ad d d V h + <f d f % C V C = g ad f^ 
and the only non-zero coefficient is 

fjj = R(l - 3mR) . 

Hence, 

V a V b d a d b = -R(l - 3mR)d R d R 

and 

(V 6 V a ) V a ^ B = -R(l - 3mR)V b V R 4> B . 
The spinor Vr^a thus satifies 

V AA ' (VriPa) = \^ AA \ V fl ] ^ = - W a ((1 - 3mi2)V fl ^ + - , 

and <5> m d a = d R . 
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The equation satisfied by higher order radial derivatives is obtained by means of the commutator 
expansion : 



fc-1 

= £v< 

k—l 



jk-p-i 
R 



p=0 
k-l 



p=0 



(42) 



It is useful to obtain the explicit expression of the action of <E> 0a and 8r on ^a by contraction. 
We simply need the Infeld-Van der Waerden symbol : 



- A A' 

<7i 



It follows 



n\ —mi 
—rh\ l\ 



1 — 3mR a 7 1 — 2>mR ~ A i 



2/10 
F { 



2F 



The conservation law for Vr^a is the following 

<7 AA> ((V^)(V^aO) = 2K (-i?(l - 3mi?) (9 R ) a (V^a)V^) 

1 — 6mi? 



+- 



{d R ) a I/JA^R^A' 



= yj^St (2R(1 - 3mR) |(V R ^) | 2 - (1 - 6mfl)^ (V K f) j 



The I? norm of on a 2-sphere is controlled uniformally on S7+ by that of 3 Vo ■ This in turn 
is controlled by the L 2 norms of ^\ and Ori^i using the second part of Dirac's equation. And 
since 



,4 



- (4>it A D6 A + M - ^o^a) 



28 



we can estimate the L 2 norm of ipo on S 2 uniformly on S7+ Q by the sum of the L 2 norms of ip\ 
and (Vfl^)i. Then we deal with the error term as we did in section 4. 

Remark A.l. Note that the error term is much simpler than it was when we used partial 
derivatives. This is to be expected but the practical upshot here is that we do not need to control 
angular derivatives in order to get estimates on the derivative transverse to This remains 

true for higher orders. It can be useful to bear this in mind if one is interested in controlling 
transverse derivatives with low angular regularity. The same is very probably true for Maxwell. 
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